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ABSTRACT 

Since M4.5-M9 dwarfs exhibit equatorial rotation velocities of order 10 knis“^ on average, 
radial velocity surveys targeting this stellar population will likely need to find methods to effec¬ 
tively remove starspot jitter. We present the first high resolution Doppler images of the M4.5 
dwarf, GJ 791.2A, and the M9 dwarf, LP 944-20. The time series spectra of both objects re¬ 
veal numerous line profile distortions over the rotation period of each star which we interpret 
as starspots. The transient distortions are modelled with spot/photosphere contrast ratios that 
correspond to model atmosphere temperature differences of Tphot - Tgpot = 300 K and 200 K. 
GJ 791.2A is a fully convective star with v sin i = 35.1 kms“^. Although we find more starspot 
structure at high latitudes, we reconstruct spots at a range of phases and latitudes with a mean 
spot filling of ^ 3%. LP 944-20 is one of the brightest known late-M dwarfs, with spectral type 
MOV and v sin i = 30.8 kms“^. Its spectral time series exhibits two dominant transient line 
distortions that are reconstructed as high latitude spots, while a mean spot filling factor of only 
1.5% is found. The occurrence of low-contrast spots at predominantly high latitudes, which we 
see in both targets here, is in general likely to be responsible for the low amplitude photometric 
variability seen in late-M dwarfs. For GJ 791.2A, the radial velocities induced by the starspot 
features yield an r.m.s. velocity variability of 138 ms“^, which can be reduced by a factor of 1.9 
using our reconstructed surface brightness distributions. 

Subject headings: stars: individual(GJ 791.2A, LP944-20) — methods: data analysis — techniques: 
spectroscopic — stars: late-type — starspots 1 


1. Introduction 


The Kepler mission has yielded many tran¬ 


siting planet cand idates ([Borucki et al. 12011 


iBatalha et alJ 12013 ) while ground based radial 
velocity (RV) surveys at r e d or infrared wave¬ 


lengths (jBean et al.l 


vs at rea 

J lioTot E 


arnes et al.l l 2012 [ ) 


have now demon strated precis i on in late-M stars 
of ^ 2.5 ms“^ ( Barnes et al ] HE! indicating 
the potential to detect rocky planets down to 
M ^ 1.5 M 0 . For early-M stars, planet occur¬ 


rence fractions of order 0.5 (iBonfils et al.l 120131) 
and even > 1 per star (|Tuomi et al.l l2014l) are 


expected, so understanding any limitations from 
stellar variability is vital. Meanwhile, M dwarf 
planet hosting stars have been found to exhibit 
significant ac tivity with photometric variability of 
up to ^ 10 % ( Knutson et al.ll 201 ll : ICrossfield et ah 
[ml. Activity will limit searches for low-mass 


planets around low mass stars 
number of epochs ([Barnes et al 


m 


a r easonable 
201l[ ) because 


by M4.5 y stars on average become moderate ro¬ 


tators ([Jenkins et al 
^ 5 kms“^ 


2009[ ). with mean 


V sm % 

while rotation is even greater 15 


kms) for the latest-M stars despite a decline in Ha 
activ ity ( Mohantv fc Basiil 2003 : Reiners fc Basril 

Iml. 

Our understanding of magnetic field generation 
in fully convective stars has only seen advances 
in relatively recent times. iReiners fc Basril ([2007h 


first observed strong magnetic fields in fully con¬ 
vective stars from Stokes I measurements of FeH 
lines. A further study ( Rniners fc Basril 120091 ) 
found that only 14% of the total magnetic flux 
detected in Stokes I is seen in Stokes V. Because 
Stokes I sees the total magnetic field and Stokes 
V only sees the net magnetic field, this implies 
that most of the field is arranged on small scales 
over which opposite polarities cancel each other 
out. The net field in early M dwarfs, which 
still possess radiative cores, is nevertheless even 
lower, at only 6 %. This implies that fully con¬ 
vective stars possess more ordered dipolar field 
structure in agreement with the Stokes V Z eema n 
Doppler imaging st udies of IPonati et al.l (120081) , 


Morin et al 


and Morin et al.l (120101) . Do- 


nati et al. and Morin et al. found that the topol¬ 
ogy amongst the earliest M stars (MOV - M3V) 
was dominantly toroidal with non-axisymmetric 
poloidal field structures, while the mid-M (M4V) 


stars at the fully convective transition bound¬ 
ary exhibited axisymmetric large-scale poloidal 
fields. The latest stars in the sample (M5 - M 8 V) 
were found to possess either strongly axisymmet¬ 
ric dipolar fields or weak fields with a significant 
non-axisymmetric component. The stars in these 
studies did not however possess sufficient rotation 
to enable brightness imaging. 

Applying Doppler ima ging techniques to rapidly 
rotati ng M 1 - M 2 dwarfs, Barnes fc Collier Cameron 
(l 200 lh and iBarnes et al ] (l2nn4 have shown that 
early-M dwarfs are more uniformly covered with 
spots than F, G & K stars. Images of two M4V 
dwarfs howev er reveal fewer spots and low spot 
filling factors (Morin et al 


nlling lactors (IMorm et al.ll^UU? 
20091 ) . While ICrossfield et al 



Phan-Bao et al.l 


20141 ) have pub¬ 
lished a cloud map of the cool brown dwarf, Luh- 
men 16B, no images of late M dwarfs have so far 
been published. In this paper we present Doppler 
images of the latest fully convective M dwarfs, us¬ 
ing spectral time series observations at red-optical 
wavelengths. 

1.1. GJ 791.2A 

GJ 791.2A is an M4.5 dwarf with a previ- 
Qusly estimated v sin i = 32.1 ± 1.7 kms“^ 
( Deifosse et al.l[l9^ ). At 8.84 pc it is a bright, 
nearby binary system with an astrometrically de¬ 
termined period of 1.4731 ± 0.0008 yrs, a maxi¬ 
mum separation of ^0.16^ ^ and primary compo - 


nent mass of 0.286 Mq ( Benedict et al. 2000[ ). 


Our spectra do not show evidence for the sec¬ 
ondary 0.126 Mq component, which is estimated 
to be A V = 3.23 faint e r than the primary com¬ 
ponent. iMontes et al.l ( 2001 ) found this young 
disk system does not satisfy Hyades Supercluster 
membership. 

1.2. LP 944-20 

Flaring radio emission h as been detected o n 
the M9 dwarf LP 944-20 ( Berger et al. 
which nevertheless exhibits lo w I/H^ /ThnU 
spite a V sin i = 30 kms“^ (iMohantv fc Basril 


2001 [) . 


de- 


20031 ). iRibasI ([20031 ) determined that LP 944-20 


may be a member of the Gastor moving group 
(d = 5 -20 pc fro m the Sun) with an im plied age 
of 320 ± 80 Myr. IPieterich et al.l ( 20141 ) report a 
parallax of 155.89 ± 1.03 mas (d = 6.41 ± 0.04 pc). 
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Table 1: System parameters and summary of observations for GJ 791.2A and LP 944-20. 




GJ 791.2A 

LP 944-20 

Spectral Type 


M4.5V 

MOV 

V sin i 

[kms 

35.1 ±0.4 

30.8 ±0.5 

Axial inclination 

[degs] 

54 ±9 

55 ±9 

Rotation Period 

[days] 

0.3088 ±0.0058 

0.162 ±0.005 

(derived) 

[R©] 

0.265 ±0.035 

0.120 ±0.014 

(estimated)^ 

[R©] 

0.279 ±0.028 

0.118 ±0.012 

Nobs 


111 

39 

Exposure times 

[s] 

180 

360 

Observation span 

[MJD] 

56904.13046 


2014 Sept 03 


- 56904.15523 


Observation span 

[MJD] 

56906.98665 

56907.22626 

2014 Sept 06 


- 56907.21841 

- 56907.40070 


^Baraffe (1998) assuming 320 Myr (LP 944-20) and 600 Myr (GJ 791.2A) ages. 


2. Observations 


We used the Ultraviolet and Visual Echelle 
Spectrograph (UVES) at the Very Large Tele¬ 
scope (VLT) to obtain time series spectra of both 
stars, with a central wavelength of 8300 A, and 
a spectral range of 6437 - 10253 A. Observations 
were made with a 0.4" slit, yielding a spectral 
resolution of ^90,000. We observed GJ 791.2A 
on 2014 September 3 & 6. Since pointing re¬ 
strictions were in place for the first half of the 
first night, only 11 spectra of GJ 791.2A were 
taken, spanning 0.595 hrs (35.7 mins), before tele¬ 
scope closure due to high winds. Both GJ 791.2A 
and LP 944-20 were observed on September 6 for 
5.562 hrs and 4.187 hrs respectively. Observations 
are summarised in Table [H 


Th e spectra were optimally extracted (jHorne 
19861 ) using the Starlink package, echomop. 


which removes all but the strongest sky lines. 
Error information based on photon statistics and 
readout noise is propagated throughout the ex¬ 
traction process. 


3. Doppler Imaging 

3.1. Least squares deconvolution 

To enable starspot induced line distortions to 
be detected for imaging purposes, we applied our 
i mplementation of least squares deconvolution 
( Barnes et al.lll99^ . Using this procedure, a high 


signal-to-noise ratio (SNR) mean line profile is de¬ 
rived for each spectrum with the effects of blending 
due to rotationally broadened lines removed. The 
deconvolution algorithn i and analysi s appl ied to 


(I2ni4 . We 

have used high SNR observations of non-rotating 
stars of the same or similar spectral type to de¬ 
rive empirical line lists for deconvolution. Eor 
GJ 791.2A, we used GJ 105B (M4.5V) as a slowly 
rotating standard star, while for LP 944-20, LP 
888-18 (M7.5V) was used. 

Obtaining an empirical line list involves remov¬ 
ing the blaze, normalising and stitching into a one 
dimensional spectrum. Normalisation of spectra 
is carried out by fitting a spline or polynomial to 
the spectrum. The minimum number of spline 
knots or polynomial degree to fit the continuum 
behaviour is adopted and generally varies with 
the length of spectrum being fitted. By reject¬ 
ing outliers at a specified level below each fit and 
iteratively re-fitting, a good approximation to the 
continuum level can generally be found. Eor stars 
such as M dwarfs, which contain so many lines that 
only a pseudo-continuum is seen, we found that it 
was necessary to obtain the maximum value of the 
spectrum in 10 A intervals and use these values for 
continuum fitting. 


Line wavelengths and depths are then identi¬ 
fied empirically, with normalised depths > 0.1. 
Lines with strong chromospheric components such 
as Ra and the infrared Na doublet and Ga ll triplet 
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are removed. Although LP888-18 has been used 
successfully to r ecover precise radi al velocities in 
M7- M9V stars (|Barnes et al.ll2Q14 ). a closer spec¬ 
tral type match for LP 944-20 would likely en¬ 
able line profiles with more optimal SNRs to be 
obtained. Line lists from model atmospheres do 
not provide a precise enough match to the ob¬ 
served spectra to enable deconvolution, although 
in ^3.31 we use models to obtain line equivalent 
width ratios and intensity ratios for Doppler imag¬ 
ing. By determining line lists empirically, we in¬ 
clude not only atomic lines, but also molecular 
lines, which are the dominant source of opacity 
in cooler atmo spheres. For instance, the Vienna 


Line Database ()Pisknnov et al.lll995l:lKnpka et al 


1999[ Rvabchikova et al. 2011) contains ^ 15800 
lines with normalised depths of 0.5-1.0 in the 
6500 A - 10250 A range for the coolest available 
temperature of T = 3500 K. The database cur¬ 
rently holds a limited number of molecular species, 
however for T = 3500 K ^ 95% of the opacities 
are due to TiO and MgH (comprising 92% and 3% 
respectively). Of the remaining 5% atomic metal 
lines, half have normalised line depths < 0.165. By 
3000 K, spectral energy distributions continue to 
be dominated by TiO and other diatomic molec- 
nlar bands, par t icularly CaH, VO, CrH an d FeH 
( Pavlenk3 20l4 IPavlenko fc Schmidt 2015 ). We 
estimate that 99% of opacities are due to molecules 
at T = 3000 K, which falls to 60% by T = 2300 K 
for atmospheres with 10-20% dust. The validity 
of using temperature sensitive molecular lines for 
Doppler imaging is investigated in detail in ^3.31 
For GJ 791.2A, the mean SNR of our in¬ 
put spectra was 57.5. With 7939 lines in the 
64851-102531 region, the deconvolved profiles 
possess a mean wavelength of A = 7781 A and a 
mean SNR of 3600; an improvement of ^ 63 in 
SNR. Similarly, for LP 944-20 the input and out¬ 
put SNRs were 18.9 and 584 respectively. With 
little flux in the bluest orders, we used the spectral 
range 84171-102531 (A = 9033 1) which con¬ 
tains 4150 lines, enabling an effective SNR im¬ 
provement of ^ 31 with deconvolution. During 
deconvolution, telluric lines were used to correct 
fo r drifts in UVES fo llowing procedures detailed 

(l20l4 . 


m 


Barnes et al, 


3.2. Image reconstruction algorithm 


We used the program DoTS ( Collier CameronI 


200lf ) to recover Doppler images of our targets. 


DoTS uses a model with two brightness levels to 
derive images from high resolution, high cadence 
time series spectra. A spot filling factor, /i, is 
derived for each pixel, i, on the model star by it¬ 
eratively fitting the time series spectra. 

Since there are generally more image pixels 
than data points, and the data possesses finite 
SNR, often with phase gaps, no unique image can 
fit the data for a given level of A regularising 
function is therefore imposed on the image. DoTS 
uses the image entropy 


5(/) = -E' 


i=l 


In— -h (1 
m 


M 1 - /i 

/i)ln-- 

1 — m 


( 1 ) 

which has the effect of minimising spurious cor¬ 
relations between image pixels. The function, 
5'(/), combines the entropy, / In /, of both the 
spot image, /, and the image (1 — /) In (1 — /) of 
the clean photosphere. The default model m is 
the value that a pixel will assume in the absence 
of any other constraint imposed by the data. Full 
de tails of the image recov ery procedure are given 


Collier Cameron ( 2QQl[ ). 


To obtain an image solution, the image pixels, 
fi are adjusted iteratively to maximise the func¬ 
tion 


Q{f) = S{f)-Xx^{f) (2) 

where A is the Lagrange multiplier. The value 
of the Lagrange multiplier is determined such that 
the image solution lies on the surface with 
M, where M is the number of measurements in 
the data set. In practice, systematics, or either 
underestimated or overestimated data uncertain¬ 
ties mean that that the ideal level of is not 
achieved. The point at which the rate of change 
of y^ with each iteration shows a marked decrease 
indicates a suitable stopping point. Continuing it¬ 
erations can lead to spot features being recovered 
that are not justified by the data. A fixed number 
of iterations is generally adopted for a given data 
set, which means that the optimal parameters will 
then provide the best fit (as determined by the 
minimum achieved y^) to the data. 


4 
















































GJ 791.2A - line equivalent width LP944-20 - line equivalent width 





Fig. 1.— Top: The deconvolved equivalent width as a function of limb-angle, /i, for BT Settl spectra 
compute with the WITA code (see ^3.31 for details). For GJ 791.2A (left), all spectra were deconvolved 
using a line list with T = 3000 K and for LP 944-20 (right), using a line list with T = 2300 K. Bottom: 
The center-to-limb local intensity profiles obtained from deconvolved synthetic spectra. The line profiles are 
derived from our standard stars, GJ105 (M4.5V) and LP888-18 (M7.5V). The center-to-limb variation of 
profile equivalent width in the upper panels is adopted. The intensities at each limb-angle, obtained from 
the synthetic spectra, are the mean intensities determined over the appropriate deconvolution wavelength 
ranges for our target stars. 
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3.3. Input models for imaging 

A study of G, R and I ph otometric variabil¬ 


ity in M dwarfs was made by iRockenfeller et al 


( 2QQ6h for late M dwarfs. While dust clouds 
were ruled out for the M9 star 2M1707+64, the 
light curves were best fit by a dusty atmosphere 
with cool spots and AT = lOOiC. iBerdvugina 
( 2QQ5[ ) has also shown that Tphot ~ ^spot declines 
with decreasing Tphot in giants, leading us to ex¬ 
pect contrasts corresponding to temperature dif¬ 
ferences of only a few hundred K in mid-late M 
dwarfs. Dopp ler images of V374 P eg (M4V) were 
obtained by ( Morin et al.l l2QQ8bh who employed 
least squares deconvolution using only atomic line 
lists and a blackbody temperature difference of 
AT = Tphot - ^spot = 400 K. 

Since Doppler imaging requires SNRs of a few 
hundred to identify starspots in time series spec¬ 
tra and derive reliable images, we combined the 
signal from as many lines as possible in our ob¬ 
served spectra, including both atomic and molec¬ 
ular lines. Since the behaviour of such a combina¬ 
tion of lines with different temperature sensitiv¬ 
ities is not clear, we investigated the validity of 
using least squares deconvolved mean line profiles 
for Doppler imaging of cool M dwarf atmospheres. 

Synthetic spectra were c omputed using the BT- 
Settl model atmospheres (I Allard et al.l l2Q12h bv 
the W ITA code, first described in iGadun fc Pavlenko 
(1997) using opacity sources listed in Pavlenko et al.l 
(120071) . Further details can also be found in 


avlenkol ( 20141 ) and IPavlenko fc Schmidtl ( 2015h . 
High resolution disk integrated model spectra 
were generated in the 2000 K to 3200 K in¬ 
terval in steps of 100 K. In addition, spec¬ 
tra were generated for limb-angles, /i = cos{0) 

= 0.1, 0.2, 0.4, 0.6, 0.8 & 1.0, where 0 is the limb- 
angle measured from disk center. The spectra 
were interpolated onto our observation wave¬ 
lengths to enable us to determine the contrast 
ratios and line equivalent width behaviour as 
a function of limb angle over the wavelengths 
used for deconvolution. For GJ 791.2A, we 
adopted Tphot = 3000 K and for LP 944-20, we 
adopted Tnhot = 23 00 K in accordance with 
,eid fc Hawlevl (|2005h (see also Kaltenegger fc Traub 
l20d 

In the case of GJ 791.2A, we carried out decon¬ 
volution on the synthetic spectra for T = 2600- 


3000 K (at 100 K intervals) over the same wave¬ 
length range as described in sa but using a 
model-derived line list for the disk integrated spec¬ 
trum with T = 3000 K. Deconvolution was per¬ 
formed for the spectra at each /i and the equivalent 
width of the deconvolved line was measured. We 
adopted the four-parame ter limb-darkening equa¬ 
tion used bv I Claret] ( 2000[ ) to fit the line equivalent 
width variability as a function of limb angle. In 
addition, we determined the mean intensity of the 
model spectra at each /i over the deconvolution 
wavelengths and again fit the same four-parameter 
equation to obtain the center-to-limb variation in 
intensity. We carried out exactly the same proce¬ 
dure for LP 944-20, using a model deconvolution 
line list for T = 2300 K to deconvolve spectra at 
all /i for T = 2000-2300 K. 

Figure [T] illustrates the behaviour of the line 
equivalent width and intensity as a function of 
limb-angle, ji. The equivalent width at disk center 
(/i = 1.0) increases monotonically as the temper¬ 
ature decreases from 3000 K to 2600 K. This be¬ 
haviour is the result of the line opacities becoming 
stronger at cooler temperatures. By 2500 K (not 
shown), the equivalent width begins to decrease 
again with a relative value of 1.28, (c.f. 1.26 for 
2700 K and 1.31 for 2600 K). For the cooler tem¬ 
perature range, with a deconvolution template of 
2300 K (Figure [U top right panel), a monotonic 
decrease in equivalent width at limb center is in¬ 
stead seen. This behaviour occurs because the 
line opacities are so numerous and strong that 
the lines are effectively weakened relative to the 
pseudo-continuum level. The equivalent width de¬ 
creases with decreasing limb-angle, except for T = 
3000 K, which shows a slight increase out to ji = 
0.5 (Figuredl top left). Nevertheless, at the cooler 
temperatures (Figure [U top right) this trend re¬ 
verses at the smallest limb angles for T = 2000 K 
- 2200 K. While decreases in line equivalent width 
may result from real changes in line depths, a mis¬ 
match between the line list and the spectrum to 
which deconvolution is applied will lead to non- 
optimal recovery of the line profile. This effect 
thus also contributes to trends of the recovered 
line equivalent width with limb angle and temper¬ 
ature. 

The relative center-to-limb variability is illus¬ 
trated in the lower panels of Figure [TJ The in¬ 
tensities at each temperature are relative to the 
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central intensity of the 3000 K spectrum (Figure 
[TJ bottom left) and the 2300 K spectrum (Fig¬ 
ure [U bottom right). The line profiles plotted 
in Figure [Hand used for image recovery are how¬ 
ever derived from the same slowly rotating stan¬ 
dards used in ^3.11 rather than th e model spec¬ 
tra, following iMorin et al.l ( 2008bh . Because the 
standard stars are observed with the same instru¬ 
mental setup, and because individual line opacities 
are not a good match, as noted in 113.li we expect 
that they better represent the local intensity pro¬ 
file shapes than those derived from the models. 
The equivalent width from the template star de¬ 
rived line profile is thus adopted for the imagining 
procedure while the relative equivalent width be¬ 
tween the photosphere and spot is obtained from 
the models (i.e. as indicated in upper panels of 
Figured]). 


For a star with a photospheric temperature of 
^phot = 3000 K, the line core is very close to 
the continuum level at Tgpot = 2700 K, being ei¬ 
ther above or below depending on the limb an¬ 
gle. Hence it is likely that spot features will only 
be distinguishable (in terms of contrast) if the 
temperature difference between photosphere and 
spot possess AT > 300 K. This corresponds to 
a contrast at disk center of 4^ot/4\ot = 0.41 for 
the adopted models with Tphot = 3000 K. With 
smaller temperature differences (AT < 300 K), 
we expect the contrast will be too small to en¬ 
able the spot induced distortions to be adequately 
fit with a spot filling factors <1.0. Similarly for 
a star with Tphot = 2300 K, a temperature differ- 
ence of AT ~ 200 K (/spot/.fphot = 0-48) is likely 
required to enable spot filling factors of 1.0. In 
section m we present results for image recovery 
using spot/temperature contrasts corresponding 
to AT = 300 K and AT = 200 K for GJ 791.2A 


and LP 944-20 respectively and additionally in¬ 
vestigate the effect of modifying the adopted con¬ 
trasts. 


3.4. Rotation periods and system param¬ 
eters 

The Super WASP archive ([Butters et al.l[2010l ]n 
contains 10166 observations of GJ 791.2A, with 
mean photometric errors of 0.24, spanning ^6.5 
years. Super WASP is however not sensitive to 


^http://www.superwasp.org 


objects, such as LP 944-20, with V > 15.0. We 
searched for periodicity in the data us ing a phase 
dispe rsion minimisation algorithm ( Lohr et "ahl 
2014h . obtaining P = 0.257025(1) d. The ampli¬ 
tude of variability, AV = 0.02 (8% of the pho¬ 


tometric error s), is comparable with t he ampli¬ 
tudes found bv iRockenfeller et~al] ( 20061 ) for other 
variable M dwarfs. However, folding subsets of 
data spanning only a few tens of days does not 
reveal obvious periodicity, suggesting higher ca¬ 
dence data with better photometric precision is 
needed. 

Adopting published parameters can lead to 

systematic and characteristic biases in the recon- _ 

strutted images ( Rice et al lIT^IUnruh fc Collier Garner^ 


19951). For example, the rotation velocity, v sin i, 


and line profile depth combination are corre¬ 
lated, and the optima l for a given inpu t model 
must be determined ( Barnes et al.lI2 QQq 1 ). Over¬ 
estimation of profile equivalent width typically 
leads to images with completely filled polar regions 
as the fitting procedure attempts to match the ob¬ 
served profiles by reducing the depth of the line 
center. This can be seen f or example in t he im ages 
of G dwarfs presented bv iBarnes et al. (1 19981). as 


compa red with the same images in iBarnes et al 


(|2QQll ) , where a lower equivalent width and v sin i 
were adopted following a revised parameter opti¬ 
misation procedure, which we use here. The main 
system parameters, recovered by minimising 
in the the multi-dimensional parameter space are 
listed in Table [H 

Uncertainties in parameters can be obtained 
from the likelihood function based on the mini¬ 
mum values obtained and the number of data 
points contributing to the image solution, as out¬ 
lined in a We have also considered the effect 
of systematic errors that could arise from adopting 
the wrong local intensity profiles and the AT val¬ 
ues and limb darkening parameters investigated 
in §4. For both GJ 791.2A and LP 944-20, us¬ 
ing the model-derived local intensity profiles leads 
to V sin i values that are 0.2 kms“^ lower than 
with the adopted standard star local intensity pro¬ 
files for instance. These additional sources of 
systematic error are included in our error esti¬ 
mates in Table 111 indicating v sin i uncertainties 
of ^ 0.5 kms~^ and inclinatio n uncertainties of or¬ 
der 10°. Morin et al. ( 2008bh similarly estimated 
V sin i uncertainties of ^1 kms“^ for the M4V 
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Fig. 2.— Reduced (Xr) plots of rotation period, P, vs axial inclination, i with 67.3% confidence level. 
Xrmin = 1-354 (GJ 791.2A) and x^min = 1.173 (LP 944-20). The derived parameters and uncertainties are 
listed in Table 1. 


star, V374 Peg. 

The axial rotation period 


P = 


27tR^ sin i 


V sin i 


(3) 


can reliably be determined by utilising re¬ 
peated starspot distortions in the line profiles 
for > 1 stellar rotation. We assumed uncertain¬ 
ties of 10% for the estimated radii of GJ 791.2A 
and LP 944-20 in Table 1, which are in agree¬ 
ment with the discrepancies between theoretical 


2012 ; 


and observed M dwa rf radii (|Stassun et al 
Williams et al. l2ni,^ . Although GJ 791.2A is 
considered a young disk object, which was not 
found to be a member of the Hyades superclus¬ 
ter by iMontes et al.l ( 200lh . we have assumed a 
radius appropriate f or a 600 Myr object following 
Baraffe et al.l (|l998h . Our radius uncertainty of 
10% c orresponds to a minim um age of 100 Myr. At 
3 Gvr. lBaraffe et al] (|l998h estimate a radius that 
is < 2% larger than at 600 Gyr, which is some¬ 
what lower than our ado pted 10% uncertain ty. 


From the estimated radii (|Baraffe et al.Nl998[) in 
Table 1, we find maximum periods for our tar¬ 
gets of Pestimated/sini = 0.4023 ± 0.0406 d and 
0.1939 ± 0.0200 d respectively. Using DoTS, we 
then performed x^ minimisation searches for P vs 
i (for P < Pestimated), which yielded P = 0.3088 
d and i = 54° for GJ 791.2A. Similarly for 
LP 944-20 P = 0.1620 d and i = 55° (Figure 


ED- That periodicities can be determined this 
way indicates coherence of spots on a single ro¬ 
tation timescale for LP 944-20 and on three day 
timescales for GJ 791.2A. We note that the Su- 
perWASP period of 0.257 d leads to a significantly 
poorer fit to the data with mis-match of repeated 
starspot features. The derived radii for both stars 
are consistent with the estimated model radii with 
the 10% uncertainties. 

4. Results 

The observed time series deconvolved profiles, 
along with the regularised fits using the param¬ 
eters in Table [U are shown for both stars in 
Figure m To enhance the appearance of starspots, 
the time series profiles are shown as 2D colorscale 
plots divided by the mean profile in Figure |4] for 
GJ 791.2A and Figure [5] for LP 944-20. Figures 
also show the fits with the mean profile 
removed along with the residual time series. The 
reconstructed maps in Figures [6] - [HI which have 
2° pixel resolution, reveal spots reconstructed at 
a range of latitudes. The gradient of a starspot 
feature in the time series spectra determines the 
stellar latitude at which it is reconstructed, while 
the recovered phase is given by the time at which 
the feature crosses the profile center. The right 
hand panels in Figures [6] - |H1 indicate the mean 
spot filling at each latitude, while the mean spot 
filling as a function of phase is shown for different 



































Profiles & Fits - GJ 791.2A (Tspo^=2700K) 


Profiles & Fits - LP 944-20 (Tspot=2100K) 






Fig. 3.— Deconvolved profiles for GJ 791.2A (left) and LP 944-20 (right). The data are plotted as points 
along with the error bars. The maximum entropy regularised profiles corresponding to the images in Figures 
[6] & [HI are shown as line plots. Observations are phased with epochs HJDO = 2456904.63601 (GJ 791.2A) 
and HJDO = 2456907.73034 (LP 944-20), which are the mid-exposures of the first observation in each time 
series. 
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GJ 791.2 - observed 



Fig. 4.— Mean profile divided deconvolved line profile time series (left), fits (middle) and fit residuals (right) 
for GJ 791.2A using a spot temperature of Tgpot = 2700i^ (upper panels) and Tgpot = 28007^ (lower panels). 
The starspot features appear as white trails while black trails are regions that are more negative relative to 
the mean profile. The spectra are centered on the observed rest frame and phased according to the adopted 
0.3088 d period of GJ 791.2A. The vertical lines indicate a v sin i profile width of of 35.1 kms“^. 
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Phase 




Fig. 5.— As for Figure IH but for LP 944-20. 
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Fig. 6.— Maximum entropy regularized image reconstructions of GJ 791. 2A plotted as Mercator (top) and 
3D projections (bottom). The image of GJ 791. 2A is shown for parameters optimised with Tphot = 3000 K 
and Tspot = 2700 K, which have a contrast ratio at disk center of /spot(0)//phot(0) = 0.41. The tick marks 
in the Mercator plot indicate the phases of observation. The latitudinal mean spot filling is shown in the 
upper left panel. The mean phase spot filling is indicated in the panel below the Mercator map for latitudes 
-30° to 90° (solid line), -30° to 45° (dotted line) and 45° to 90° (dashed line). 


latitude ranges in the panels below each Merca¬ 
tor projection. The 3D images at rotation phases 
0.00, 0.25, 0.50 & 0.75 are also shown for all cases. 

4.1. GJ 791.2A 

GJ 791.2A shows a number of starspot features 
of different sizes distributed across most phases. 
We used Tphot = 3000 K and Tgpot = 2700 K, with 
a disk center contrast ratio of Jg^pot/^phot ~ 
for the image reconstruction in Figure [6l The 
phase range (j) = 0.0000- 0.0803 is provided by 
observations on September 3 that align spot fea¬ 
tures with those from September 6, which first 
come into view at phases 0 ^ 0.8-1.0. Phase over¬ 
lap between nights is confined to the narrow (p = 
0.0000-0.0008 range, corresponding to the first 
observation on September 3 and the last observa¬ 
tion on September 6. Any starspots at the center 
of the 0.0803 - 0.2504 phase gap in observations 
are still seen for 66% of the total time they are 
physically visible on the star. Spots in the phase 
gap can thus be reconstructed, albeit with poten¬ 
tially decreased reliability in latitude. 


Since cos(6>) oc starspot trail gradient to first or¬ 
der (where 0 is the stellar latitude), and cos(6>) is 
a slowly changing function when 0 is small, the 
elongation of spots at low latitudes is a refiec- 
tion of the uncertainty in 0. There are ^ 17.5 
resolution elements in the rotationally broadened 
profile of GJ 791.2A, equivalent to ^ 10.3° res¬ 
olution, or 0.029 in phase, at the equator (the 
180 sec exposures result in rotational blurring of 
only 2.2°). The weakest low latitude features ap¬ 
pear elongated with FWHM values of A0 ^ 4° - 6° 
and A6> ~ 22°-30°. These features are approx¬ 
imately equivalent to circular spots with radii 
- 9.5 - 13.5°. 

The maximum spot filling factor in the image of 
GJ 791.2A is /max = 0.823 (82.3%) and the mean 
spot filling is / = 0.032 (3.2%). Spot activity is 
seen at all latitudes, but with peak filling factors 
concentrated at ^ = 59° and particularly 0 = 75° 
as indicated in Figure [6] (panel to right of Merca¬ 
tor map). Spot coverage as a function of phase 
(Figure [6l panel below Mercator map) is higher 
on average in the 45°-90° region (dashed line) 
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compared with the lower -30° - 45° region (dotted 
line), again demonstrating the greater degree of 
spot filling at higher latitudes. To some extent 
this is due to the fact that the number of reso¬ 
lution elements is greater at lower latitudes where 
features are better resolved. In other words, in the 
-30° - 45° range, the spot coverage as a function of 
phase approaches zero in regions outside the spots, 
which can be seen to produce more narrow peaks 
than for the 45° - 90° region. 

4.1.1. Images with AT = 200 K & 4OO K 

We also attempted to fit our spectra with 
^spot = 2800 K (AT = 200 K) with an inten¬ 
sity contrast ratio at disk center of 
= 0.57. However, for ^spot = 2800 K, we achieved 
a reduced of Xr — 2.39, whereas with ^spot 
= 2700 K, we achieved Xr — 1-35. The contrast 
in the corresponding fits is lower in Figure |4] (lower 
panels) with noticeably greater residuals. The up¬ 
per panels in Figure [3 show the image for T^pot 
= 2800 K which shows a maximum spot filling of 
/max = 0.983 and a mean spot filling, / = 0.038. 
The image is very similar to the Tgpot = 2700 K 
image, except for the greater degree of spot filling. 
As expected, DoTS has tried to account for the 
lack of contrast by enlarging some of the spots 
(e.g. the spot at phase 0.62 and -30°) in an at¬ 
tempt to fit the starspot distortions in the line 
profiles. 

With Tgpot = 2600 K (AT = 400 K), we find 
a contrast of I^^ot/^phot ~ 3?he recon¬ 

structed image closely resembles the AT = 300 K 
image, but shows reduced spot filling factors of 
/max = 0.690 and a mean spot filling, / = 0.026 
(x? = 1.38). 

4.2. The effect of limb darkening uncer¬ 
tainty on image reconstruction 

The concentration of spots seen on GJ 791.2A 
in Figure [6] at 75° appear to occur in a relatively 
narrow band, which is most evident in the 3D plots 
shown below the Mercator projection in Figure 
[6l Features that do not traverse the whole v sin i 
extent of the star are reconstructed at high lati¬ 
tudes. Closely separated starspot trails (in phase), 
or broader, more complex trails can not be re¬ 
solved and can result in merging of starspot fea¬ 
tures in the images. At higher latitudes, this can 


lead to the appearance of arcs or rings of structure. 
To further investigate the axisymmetric structure 
at 75° on GJ 791.2A further, we have explored a 
number of different input models. 

Overestimating the degree of limb darkening 
leads to a more V-shaped line profile, requiring a 
higher v sin i to fit the time series spectra. On the 
other hand, a lower degree of limb darkening leads 
to a more U-shaped profile and a typically lower 
V sin i estimate, potentially reducing circumpolar 
structure. To assess the reality of the high latitude 
axisymmetric structure on GJ 791.2A, we investi¬ 
gated the effect of modifying the limb darkening 
in our input models. We scaled the center-to-limb 
intensities at each limb-angle using a linear rela- 
tionship, such that /adjusted = /model(l-^(1-M))> 
where ^ is the scale factor. We investigated ^ = 
-0.2, -0.1, 0.1 & 0.2 and found that the x^ 
fitting our profiles was minimised by reducing the 
limb darkening by 10% (i.e. ^ = —0.1). A re¬ 
duced Xr — 1*13 was achieved compared with Xr 
= 1.35 for the default case. The resulting image is 
shown in Figure El (lower panels) which shows that 
the axisymmetric ring at high latitude structure is 
less prominent. Mean and maximum spot filling 
factors of / = 0.033 and /max = 0.76 are found. 

With a 20% reduction in limb darkening, the 
high latitude axisymmetric ring structure is re¬ 
duced in prominence, leaving only the main con¬ 
centrations of spots, but all lower latitude struc¬ 
ture becomes confined to the latitude range 0° - 
30°(x^ = 1.38 is achieved). Increasing the degree 
of limb darkening led to strong polar filling factors 
and poorer fits. We also investigated using Tphot 
= 3100 K with AT = 400 K (a contrast ratio of 
^spot/^phot ~ 0.45). In this case, the center-to- 
limb intensity variation was similar to our opti¬ 
mal limb darkening (i.e. when reduced by 10%) 
for Tphot = 3000 K. The center-to-limb trend is 
however more linear and leads to an augmenta¬ 
tion of the axisymmetric ring compared with the 
image in Figured with a reduced Xr = 2.04. 

Finally, removing the center-to-limb line equiv¬ 
alent width variation had little discernible effect 
on the image, although the strength of the spot 
filling at 75° increased by 9%. Optimising the 
center-to-limb variation in intensity and the non¬ 
linearity of the limb darkening is thus important 
for ensuring that artefacts are minimised. The 
need to reduce the limb darkening by 10% gives 


13 


an empirical indication of the systematic uncer¬ 
tainty in the current models for cool M dwarfs. 

4.3. LP 944-20 


as for our default AT = 200 K, but with en¬ 
hanced circum-polar filling, as might be expected 
by increasing the limb darkening (/ = 0.019, 

/max = 0.88). 


The two dominant starspot trails seen in the 
time series spectra of LP 944-20 (Figure [5]) pos¬ 
sess steep gradients, leading us to expect starspots 
should be reconstructed at predominantly high 
latitudes. The image of LP 944-20 in Figure [8] 
(top two panels) is reconstructed with i = 55° and 
P = 0.1620 d and shows activity is concentrated in 
two groups at 0 ^ 0.2 & 0.65, approximately half 
a rotation phase apart. For v sin i = 30.8 kms“^, 
the equatorial spot resolution is 11.5°. As with 
GJ 791.2A, spots are found at high latitude, with 
a maximum in the azimuthally averaged spot fill¬ 
ing at 79°. The spot reconstructed at 0 = 0.16 
and 0 = 55° is predominantly responsible for the 
mean latitude filling peak at ^ = 55°. Given the 
lower S/N of the data, the smaller, weaker spots at 
(j) = 0.65, 0 = 55° and 0 = 0.925, 0 = 35° may be 
artefacts in the reconstruction. We also note that 
if lower latitude spots of similar size to those seen 
on GJ 791.2A are present on LP 944-20, they may 
not be recovered owing to the lower SNR of the 
data. LP 944-20 reveals a very low mean spot fill¬ 
ing factor of / = 0.015 (1.5%) for the image recon¬ 
structed with Tphot = 2300 K and Tgpot = 2100 K 
(AT = 200 K), corresponding to a contrast ratio 
at disk center, ^spot/^phot ~ 0.48. A maximum 
spot filling factor, /max = 0.766 (76.6%), is found 
and we achieve = 1.17. 

4 . 3 . 1 . LP 944.-20 image with AT = 100 K 

The image reconstruction with Tgpot = 2200 K 
(AT = 100 K, /^potZ-fphot = 0-64) is shown in Fig¬ 
ure |8] (bottom two panels). As with GJ 791.2A, 
when we investigated a reduced contrast, we find 
that the data are not fit so well (y^ = 1.22) com¬ 
pared with Tgpot = 2100 K, and that features are 
reconstructed with larger areas to account for the 
lower contrast ratio. As expected, the mean and 
maximum spot filling values are higher than for 
(AT = 200K), with / = 0.017 and /max = 0.866. 
The lower SNR of the LP 944-20 time series makes 
investigation of limb darkening modification more 
difficult. However, we found that even for ^ = 0.2 
(20% increase in limb darkening), a = 1.17 
is achieved. The image shows the same features 


5. RV variability 


We have cross-correlated the observed time se¬ 
ries to determine the RV variability induced by the 
spots on both targets. We find a semi-amplitude 
of = 160 ms“^ for GJ 791.2A, with a corre¬ 
sponding RV r.m.s. of 138 ms“^. We also cross- 
correlated the time series line profile fits for AT 
= 300 K. By subtracting the RVs derived from 
the fits from those mea sured for the data, follow¬ 
ing IPonati et al.l (|2014| ) , the starspot induced ra¬ 
dial velocity variability in the data can be partially 
corrected. We find that the radial velocity r.m.s. 
can be reduced by a factor of 1.75, from 138 ms“^ 
to 79 ms“^. Using the optimised fits with AT = 
300 K and the limb darkening reduced by 10%, we 
find the RV r.m.s. is reduced by a factor of 1.9, to 
73 ms“^. For LP944-20, the SNR is not sufficient 
to discern regular RV variability. The radial ve¬ 
locity r.m.s. of 194 ms“^ is thus likely to be noise 
dominated. The regularised fits however yield an 
RV r.m.s value of 93 ms“^. 


Barnes et al 


found similar RV r.m.s. 
here for GJ 791.2A and 


(l2niih 

values to those found 
LP 944-20 for their scaled solar starspot model 
with V sin i = 20-50 kms“^. However, mod¬ 
els with large spot filling fractions of either 
/ = 0.30 and greater spot-to-photosphere con¬ 
trast than found here (up to Tphot/Fgpot = 0.65) 
or models with / = 0.62 and similar contrast 
(AT = 250K) are required to obtain r.m.s. RVs of 
order 100 ms“^. Based on the evidence presented 
here, it thus appears spots or spot groups are 
larger or more locali sed than the random scaled 
solar spots models of iBarnes et al.l (|201l[ ). 


6. Discussion 


The time series spectra and Doppler images 
of the two mid-late M dwarfs studied here re¬ 
veal significant spot structure, even at spectral 
type M9V. The results are in general agreement 
with the finding that t he fall-off in fractional Ha 
flux in late M dwarfs ( Mohantv fc Basrl 2003) is 


not in fact accompanied b; 
netic flux ( Reiners fc Basril 


a dr op-off in mag- 
and provide 


w a aro 

lli^, 
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Fig. 7.— As in Figure [6] for Tgpot = 2800 K (top), Tgpot = 2700 K (bottom) with limb darkening reduced by 
10% with a linear centre-to-limb scaling (bottom). 
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Fig. 8. — Maximum entropy regularized image reconstructions of LP 944-20 with panels as 
indicated in Figure [G], Images are shown for Tphot = 2300 K, with Tgpot = 2100 K (upper panels) 
and Tgpot = 2200 K (lower panels), with respective contrast ratios of /spot(0)//phot(0) = 0.48 
0.64. 
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further evidence for magnetic activity in the low¬ 
est mass M dwarfs. The image of GJ 791.2A 
shows starspots located at a range of latitudes 
and longitudes, but preferentially at mid to high 
latitudes, whereas activity is confined solely to 
high latitudes on LP 944-20. Late-F and early-G 
dwarfs, show a still greater degree of spot fill¬ 
ing at high latitudes compared with lower lati¬ 
tudes, with activity confined to localised regions 
in phase at lower lati tude (e.g. iBarnes et al. I l 2 nni[ 


Marsden et al.l I2QQ6I) . In contrast, Doppler im- 


M4V stars, the spot/photosphere contrast ratios 
are sufficient to preclude the detection of Stokes V 
signatures in the spots, meaning that magnetic 
and brightness maps indeed probe magnetic ac¬ 
tivity on different scales. 

Our images show a number of polar or circum¬ 
polar spots, especially in GJ 791.2A that might 
reasonably be expected to possess a mixture of 
magnetic polarities. This is interesting because 
the degree to which the polarities cancel will likely 
determine the large scale global field structure. 


^ _ Olectio, XXXV 

(Barnes &: Collier Cameron 2001: 

iBarnes et al. 

wnicn m luny convecrive 
erally axisymmetric and 

mia-ivi 

poloidal 

stars IS gen- 

(Morin et al. 

2004) showed distributed activitv. 

while differ- 

l2008allbl: IPhan-Bao et al. 

2009|). 

Although we 


ences are also seen when co mpared with images 


of the M4V stars V394 Peg ( 

iMorin et al. 2008b) 

and G164-31 ( 

Phan-Bao et al. 

2 OO 9 I) which possess 


similar v sin i values (36.5 kms“^ and 41 kms“^ 
respectively). Morin et al. found weak features 
concentrated at mid latitudes with only 2% mean 
spot occupancy, similar to, but slightly lower 
than the 3.3% we find for the optimised image 
of GJ 791.2A. Phan-Bao et al. on the other hand 
found only a weak de-centered polar spot with 
no low-latitude features on G164-31. The cur¬ 
rent evidence from Doppler images thus suggests 
that the lower photometric variability in cooler M 
dwarfs is the result of both lower spot contrast 
and more distributed activity at predominantly 
higher latitudes. 

Zeeman Doppler imaging does not reveal the 
field strength in the largest spots, which show sig¬ 
nificant spot/photosphere contrast, because the 
amplitude of the Stokes V signature is generally 
small. Comparing magnetic maps and brightness 
maps is often difficult as a result, and because 
Stokes V only sees the net magnetic field, thus 
making it sensitive to field structure on larger 
scales. Large scale poloidal fields were inferred 


for V394 Peg (1 a 

min et al. 2008bl) and G164-31 

Phan-Bao et al. ( 

2 OO 9 I). in ameement with the re- 


suits for more slowly rotating fully convective stars 


( Morin et al. 2QQ8al ). Both these stars showed 


only positive magnetic polarity in one hemisphere. 
These studies show contradictory results when the 
magnetic and brightness images are compared, 
since V394 Peg possessed multiple low to inter¬ 
mediate latitude features, and no polar features, 
while G164-31 showed only a single weak circum¬ 
polar feature. This probably suggests that even for 


have already noted that spots recovered in bright¬ 
ness images do not generally correlate with mag¬ 
netic field structure recovered from Stokes V, 
the lower contrasts in GJ 791.2A and LP 944-20, 
with only Tphot " Tspot = 200-300 K 
^ 0.4 — 0.5), may enable at least some of the field 
in the spots to be more readily seen in Stokes V. 
This will also depend on the mixing and degree of 
polarity cancelling on relatively small scales and 
the resolution limitations of any observations. 

With low contrast spots, future Stokes V ob¬ 
servations of M dwarfs may provide the means to 
make more direct comparisons between spot and 
magnetic field morphologies and to determine why 
the fraction of magnetic flux seen in Stokes V com¬ 
pared with Sto kes I is greater than for early M 
spectral types ( Reiners fc Basril 20091 ). In partic¬ 
ular, infrared polarimetric observations may offer 
the potential to observe the magnetic field inside 
starspots since the contrast ratio between spots 
and photosphere i s smaller at longer wa velengths. 
We also note that iHallinan et ^ ( 2Q15h have sug¬ 
gested a transition from coronal activity to auroral 
activity in late M stars and have suggested that 
such activity arises from circumpolar latitudes on 
the M8.5V, LSR J1835. It is not clear whether the 
high latitude, low contrast cool spots recovered in 
LP 944-20 could be related to, or contradict this 
mechanism. 

Dedicated RV surve ys such as the Habitabl e 
Zone Planet Finder (Mahadevan et all I2Q1Q[) . 
CARMENES (iOuirrenbach et al.ll2QlQh and SPIRou 


(|Thibault et al 

planets around 
direct influence 


2Q12h will search for low mass 
M dwarfs. Consequently, the 


of magnetic activity on stellar 
lines is of significant importance. Eor stars rotat- 
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ing with V sin i = 20 kms“^, iBarnes et al.l ( 2Qllh 
found that 5 M 0 to 20 M 0 habitable zone plan¬ 
ets could be detected around 0.1 M©, MOV stars, 
but would require 300 and 40 epochs respec¬ 
tively. Planets with 1 M 0 could nevertheless 
be detected with ^ 100 epochs for stars with 
V sin i = 5 kms“^ and (Tphot — ^spot = 250 K). 
These estimates will likely be modified by the 
findings presented here, since we find spot groups 
are more localised and larger (especially in the 
case of LP 944-20) than the spots modeled by 


Barnes et al.l (|201lf ). 


iMartfn et al.l ( 2006 ) have searched for RV vari¬ 
ability in LP944-20. In the optical, 15 observa¬ 
tions spanning 841 days, at similar instrumental 
setup to the observations presented herein, yielded 
a semi-amplitude variability of 3.5 kms“^. By 
contrast, Keck/NIRSPEC observations spanning 
only 6 nights yielded a dispersion of 360 ms“^; an 
order of magnitude smaller than at optical wave¬ 
lengths. Our finding of 194 ms“^ RV r.m.s. is 
thus closer in magnitude to the NIRSPEC ob¬ 
servations. Significant spot-coverage evolution, as 
evidenc ed through other activity in dicators in M 
dwarfs ( Gomes da Silva et al.l 2011), is thu s likel y 
to have biased the results of Martin et al.l ( 20061 ). 
Given the observed activity in M dwarfs, it is clear 
that searches for low mass planets will benefit from 
techniques that can effectively reduce or remove 
the substantial RV noise. Although one might not 
routinely expect to observe stars with rotation ve¬ 
locities of order 30 kms“^, the fact that Doppler 
imaging can be used to mitigate the effects of spots 
on GJ 791.2A by a factor of ^ 2, indicates that it 
may also be a useful method for slower rotators. 
This is particularly important in mid-late M stars, 
since they show equato rial rotation velocit ies of 
^10 kms“^ on average (| Jenkins et al. 2009h . Ob¬ 
serving only the very slowest rotators is likely to 
bias the sample to low axial inclinations, for which 
it would be much more difficult to detect orbiting 
planets. 


7. Summary 


trast at disk center, = 0.41) is required 

to model spots on the M4.5V, GJ 791.2A, whereas 
the data can be fit with AT = lOOK- 200 K (con¬ 
trast at disk center, ^spot/^phot ~ 0.64-0.48 for 
the MOV LP 944-20, albeit with a marginally bet¬ 
ter fit for AT = 200K. The factor of two lower 
spot filling for LP 944-20 (/ = 0.015) compared 
with GJ 791.2A (/ = 0.033) suggests a decline in 
magnetic activity by MOV, as seen from other 
magn etic activity indicators ( Mohantv fc Bas^ 
I2OO3I ). However, Doppler imaging at optical 
wavelengths is extremely challenging for even 
the brightest MOV targets and it is possible that 
smaller, low latitude features have not been recov¬ 
ered on LP 944-20, owing to the lower SNR of the 
data. Our ability to reconstruct small scale bright¬ 
ness features at all latitudes should benefit from 
the next generation of i nstruments cofferin g higher 
throughput and optical ( Zerbi et~al]l2014{ ) to near 
infrared wavelength coverage in conjunction with 
near infrared polarimetric measurements.. 
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